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Abstract 

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for 
| semilinear parabolic initial boundary value problems. Unlike any classical counterpart, this class 

l offers arbitrarily high, optimal order convergence. In support of this claim, error estimates are 

proved, and computational results are presented. Furthermore, it is noted that special Runge- 
Kutta methods allow computations to be performed in parallel so that the final execution time can 
be reduced to that of a low order method. 
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1 Introduction 

In this paper, three classes of semilinear initial boundary value problems are considered. Specifi- 
cally, the goal is to construct and analyze fully discrete approximations to the unique solution u(x, t ) 

dtu = —Lu + f in fix [0,(*] 

i u = 0 on an x [0,t*] 

u(x,0) = u°(x) in n, 

N 

Lu = - Y, a if (4y(x)a xy u) + £o(x)u, 

«»J = 1 

and / takes one of the following forms: 


of: 

(i.i) 

where: 


(l.l.i) f = f(x,t,u), 

(l.l.ii) / = g(x,t,u)- Vu, 

(l.l.iii) f = f(x,t,u,Vu). 

Here, n is a bounded domain in H N with dfl sufficiently smooth. Also, Ay(x) and £o(x) are assumed 
to be smooth. Further, on fi, the matrix {Ay} is symmetric and uniformly positive definite and 
Iq is nonegative. The initial data u° is assumed to be both sufficiently smooth and compatible, 
and precise hypotheses on the required smoothness of the solution u are made as needed. Then, 
in the respective sections, it is assumed that there exist constants p, c p > 0 such that for all 
(x, t) G n X [0,t*], one of the following local Lipschitz properties holds. Specifically, it is assumed 
that in case (l.l.i): 


(1.2) \u{x,t)-U\<p =>• \f{x,t,u(x,t))~ f{x,t,U)\<c p \u{x,t)-U\, 

in case (l.l.ii), for 1 < t < N: 


(1.3) \u{x,t)-U\<p =>► \gi(x,t,u(x,t))- gi{x,t,U)\<c p \u(x,t)-U\, 

and in case (l.l.iii), for 1 < j < N: 


(1.4) 


|u(x, t ) - U 0 \ < p 
and 


max |d z .u(x, t) - UA < p 
l<i<N 1 


|/(x,t,u(x,f),u) - f(x.,t,Uo, U)| < c>(x,f) - Uo\ 

and 

| f(x,t,U 0 ,V j ) - f{x,t,U 0 ,V)\ < c p \d x .u(x, t) - Uj\, 


where Vf = Ui,i ± j and Vj = d Xj u[x, t). In addition, to prove certain estimates for case (l.l.iii), 
it is assumed that: 


(1.5) 


|u(x,<) -U\<p 
and 

max |d z .u(x,t) - UA < p 
l<i<N' * v ' 


max |a^ dy /(x,t,[/o, u )| < c p 


1 <i,i<N 
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where d^f denotes partial differentiation of / with respect to its argument connected with the jth 
spatial derivative of tt. See the remarks prior to Proposition 4.3 for an explanation of the division 

(l-l-i) - (1-1-iii). 

Now, for 1 < p < oo and integers s > 0, let W* tP = W^ s,p (f2) represent the well-known Sobolev 
spaces consisting of functions with (distributional) derivatives of order < s in L p = L p (f2). Also, 
let || • || w*p denote the usual norm. Then, in particular, take H 8 = W 8 ' 2 and denote its norm by 
|| • || 5 . In addition, let Hq be the subspace of H l consisting of functions vanishing on dQ in the 
sense of trace. Its dual is denoted by H~ x with norm || • ||_ 1# Next, let the inner product on Li be 
denoted by (■,•), and the associated norm by || • ||. Further, || * represents the norm on L 0 0 , 
and || • || a>00 the norm on ioo([0, t*], H 8 ). See Adams [l] for more details. 

Equipped with the above notation, let the following be combined with (1.3) and (1.4), respec- 
tively. Specifically, assume that in case (l.l.ii): 

(1.6) max sup ||ffi(*,u(t))|| w i,oo = c g < oo, 

l<*<AT0<t<t* 

and in case (l.l.iii): 

(1.7) max sup \\d d .f(t, u(t), Vu(i))||^i,«, = c f < oo. 

l<J<N 0 <t<f 

Now, let L be extended to have domain H 2 CiHq. Then, L is i 2 -selfadjoint and for every nonegative 
integer s, it is bounded from H 8+ 2 n Hq into H 8 . Furthermore, introducing the solution operator 
T for the elliptic problem: 

{ Lv = w in fi 
v = 0 on 9fl 

as Tw = v, it is well-known (Friedman [12]) that for every nonnegative integer s, T is bounded 
from H 8 into H ,+2 n Hq. Also, the solution operator is positive definite and selfadjoint on L<l\ 
hence, T has a square root and it can be shown (Thomee [18]) that: 

(1-8) \\Th\\ < c 1 1 v 1 1 — ! Vv€ H~ l . 

Note that here and throughout this work, c (sometimes with a subscript) is used to denote a general 
positive constant, not necessarily the same in any two places. Moreover, if in a given (in) equality, 
there is a crucial element upon which c is meant to depend, such dependence is indicated explicitly. 

A rough description of the results now follows. For this, let h and k denote spatial and temporal 
discretization parameters respectively, and suppose that U£ is a fully discrete approximation to 
u(nfc) obtained according to (1.34) described below. Now, in section 2, the error committed for the 
approximation of the solution to (1.1) in case (l.l.i), is shown to be of optimal order in L 2 : 

(1.9) max || U£ - u"|| = 0{h r + k u ) 

under the condition that h~ N ! 2 {h r + k v ) is sufficiently small. Here, r and v represent respectively, 
optimal exponents, characteristic of the Galerkin method and the Runge-Kutta method upon which 
the fully discrete scheme is based. Next, section 3 deals with case (l.l.ii) and the same optimal 
Li estimate is established but under the additional condition that h~^l 2 k~ x l 2 {h T + k v ) is small 
enough. Finally, in section 4, case (l.l.iii) is studied, and it is proved that the error is of optimal 
order in H l : 
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(1.10) max || US - u n ||! = 0 {h'' 1 + k u ) 

provided h~ N l 2 {} i r_1 + k u ) is small enough. Then, a duality argument is used to obtain (1.9). For 
each of the cases (1.1. i) - (l.l.iii), results for the starting scheme (1.37) are stated without proof, 
since what is presented for the principal scheme (1.34) captures the main ideas with fewer details. 
Nevertheless, complete proofs are provided in [14]. Also, in the latter, linear problems with time 
dependent coefficients and quasilinear problems are considered. Further, preconditioned iterative 
methods are combined with specially constructed Galerkin/Runge-Kutta schemes and the results 
obtained are similar to those reported here. 

It should also be mentioned that the discovery of the methods described below was fortuitous. 
Note that there are extrapolation options other than (1.33) and (1.36), which are apparently more 
natural. For example, T*e could be used instead of /! A l e in (1.32) and (1.35), since this is suggested 
by the case that / = /(x, t). On the other hand, an iterative procedure could be used during 
initial time steps to approximate fully implicit stages, and this might be followed with a standard 
extrapolation using previously computed stages to approximate stages for the current time step. 
Both of these ideas are considered in a computational section. However, under rather general 
conditions, optimal order convergence is proved and demonstrated computationally only for (1.34) 
and (1.37). In fact, for the linear nonhomogenous problem, Crouzeix [8] has developed an explicit 
example showing that unless d|u is in the domain of L m for certain / and m, a classical fully discrete 
scheme fashioned after (1.25) cannot be expected to offer optimal order convergence. 

In [4], Baker, Dougalis, and Karakashian analyze Galerkin/Multistep fully discrete approxi- 
mations for the solution of (1.1) in case (1.1. i). Assuming local Lipschitz properties, they prove 
an optimal L 2 error estimate such as (1.9). In connection with case (l.l.ii), they have obtained 
similar results for the Navier-Stokes equations [3]. Also, for quasilinear equations which embrace 
case (l.l.ii), Bramble and Sammon* have announced results including an optimal L 2 estimate for a 
Galerkin/Obrechkoff method which is fourth order in time. Finally, note that in [11], Dougalis and 
Karakashian analyze Galerkin/Runge-Kutta approximations for the Korteweg-De Vries equation, 
and they prove optimal L 2 estimates for some modified IRKM’s which are up to fourth order in 
time. Hence, the spirit of their work is similar to that of the present study. 

In the remainder of this section, there is a presentation of material relevant to the spatial and 
temporal discretizations considered here, which concludes with a precise definition of the schemes 
for which the above claims are made. 


Spatial Discretizations 

In terms of the solution operator, (1.1) can be written as: 


(1.11) 


d t Tu = -u + Tf 
u(0) = u°. 


For the spatial approximation of the solution to this problem, let {Sh}o<h<i be a family of finite- 
dimensional subspaces of H 1 . Then, suppose that a corresponding family of operators {TX}o<h<i 
is given satisfying: 


‘Bramble, J. H., and Sammon, P. H., “Efficient Higher Order Single Step Methods for Parabolic Problems: Part 
II.” 
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i. Th : L 2 — ► Sh is selfadjoint, positive semidefinite on L 2 , and positive definite on S/,. 

ii. There is an integer r > 2 such that: 


(1.12) ||(T — Th)v\\ + h||(T — Tfc)v||i < ch 4 ||v|| s _2 Vt >eH- 2 , 2 < s < r. 

Now, problem (1.11) has the following semidiscrete formulation. Find Uh : [0, t*] — ► Sh such that: 

,o\ / dtThUh = -Uh + Thfh 

(113) 1 «*( 0 ) = ul 

where fh represents / depending on instead of u as indicated in one of (1.1. i) - (l.l.iii), and 
G Sh is a suitable approximation to u°. 

To make the machinery more definite, consider the following Ordinary Galerkin Method. From 
(1.1), let £>(•,*) be a symmetric bilinear form defined by: 

N 

D ( v > w ) = X) d x .w) + (£ov, w) V,we Hq. 

t,i=l 

Then, take Sh to consist of continuous, piecewise polynomials of degree < r — 1, vanishing on dtt. 
Now, let Th : — * Sh be defined by: 

D{T h w,x ) = (t v,x) Vtn G L 2 , Vx € S h . 

For more examples of Galerkin methods satisfying the assumptions above as well as others below, 
see Bramble, Schatz, Thomee, and Wahlbin [5], and the references cited therein. 

Next, the following inverse properties are prescribed for 5^. Throughout this work, it is assumed 
that Sh C Loo and: 

(1-14) ||xiu„<c/r"/ 2 ||x|| v X €5 fc . 

Moreover, for cases (1.1. ii) and (l.l.iii), it is assumed that Sh C W 1,<x> and hence: 

(115) Hx||wm.°° < cJT^/ 2 Hxlli VxG5 fc . 

In fact, for certain estimates related to case (l.l.iii), it is assumed that: 

(116) ||xlki.F <^-^*-^11x111 Vx€S„, p> 2. 

For details connected with (1.14) - (1.16), see Ciarlet [7]. 

According to the properties prescribed above, the restriction of Th to Sh is invertible and its 
inverse is henceforth denoted by Lh . Since Lh is also positive definite and selfadjoint on Sh, both 
Lh and Th have square roots but it is also assumed that: 

(117) IIT^IU < c|H| Vt veL 2 , 

and: 
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( 1 . 18 ) llxlli < c||L| x || v x es k . 

Defining the elliptic projection operator as Pg = T^L, it follows from (1.12) that: 


(1.19) ||(/ - PjsMI + A||(7 - 7fe)t>||i < eh*\\v[ 


Vv€H a nH%, 2<s<r. 


Also, it can be shown that L^Th = Pq is the orthogonal projection of onto S^. Then, since 
7 — Po is majorized by / — Pg in Li, it follows from (1.19) that: 

(1.20) ||(/ - PoH| < ch*||v|| 4 Vt> G H s D Hi, 2 < s < r. 

In addition to (1.19), with oj(t) = Pgu(t), it is assumed for case (l.l.i) that: 


( 1 . 21 ) 


sup yu(t) - w(t)!! ioo = 70(^1) -*■ 0, 

0<t<t* 


Further, for cases (1.1. ii) and (l.l.iii), it is assumed that: 


( 1 . 22 ) 


sup ||u(t) - w(t) 11^1,00 = 71(h) -> 0 , 
0 <t<r 


as h — * 0. 


as h — ► 0. 


For details connected with (1*21) and (1.22), see Rannacher and Scott [16], and Schatz and Wahlbin 

[17]- 

Now, problem (1.13) takes the following form. Find : [0 ,f*] — > Sh satisfying: 


(1.23) 


f d t Uh = -Lhtih + P 0 f h 
\ u fc (0) = u° h . 


In [4], Baker, Dougalis and Karakashian analyze approximations of the form (1.23) with / as in 
(l.l.i). Assuming local Lipschitz properties, they prove an optimal Li estimate: 

||« - «fc||o,co = 0(h r ). 

In connection with case (l.l.ii), Baker* has obtained similar results for the Navier-Stokes equations. 
Then in [10], Dendy studies Galerkin semidiscrete approximations for quasilinear problems which 
embrace case (l.l.iii). Assuming global Lipschitz properties, he establishes the following: 


|« - Ufclkoo = 0(h r *) 


|«-U h ||o,oo = 0(h min ( r - 2r - 2 - Ar / 2 )). 


In the present paper, semidiscrete approximations are not analyzed. Instead, (1.23) serves only as a 
source of inspiration for fully discrete approximations, and is not even mentioned in forthcoming 
proofs. 

Temporal Discretizations 

For the temporal approximation of the solution to (1.23), Implicit Runge-Kutta Methods 
(IRKM’s) are now introduced. Given an integer q > 1, a q-stage IRKM is characterized by a 
set of constants: 

* Baker, G. A., “Galerkin Approximations to the Navier-Stokes equation.* 
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“1« 


T 1 


a 9 i 

til 


a g<? 

fei. 


and it is convenient to make the following definitions: 


A = { a ij}i<i,j< q , T = diag {r,}, b T = <fei,fe 2 , . . . ,6,), e T = (1, 1, . . . , 1). 

1<I<9 

For the IRKM formulation used in this work, choose arbitrarily, to G R, yo € R", F : R n+1 — > R n 
sufficiently smooth, and k > 0 sufficiently small, so that for to < t < to + k, smooth functions 
y,y : R — > R" are well-defined by: 


(1.24) 

and: 


(1.25) 


( D t y(t) = F(t,y(t)) 

\ y(*o) = yo, 

y’W = yo + (*-*o)X! a *'i F ( t o + T- J (t-to),y y (t)), l < i < ? 

9 

y(t) = yo + (t - to)J^biF(t 0 + Ti(t - to),y*(f))- 


The method is described as explicit if a,j = 0, t < j and implicit if for any t , a,-,- ^ 0. Also, it is said 
to have order v if for every y and y defined as above, Z?fy(to) = -D{y(*o)> 0 < / < v. Butcher [6] 
has developed simple conditions for the above parameters, which guarantee a given order; however, 
only the following is explicitly required in this work: 


(1.26) l\b T A l ~ l e = 1 1 </<i/. 

To see the roots of condition (1.26), let (1.24) have n = 1, to = 0, yo = 1, and F(y) — -y, so that 
y(t) = e -t . Then, from (1.25), y(t) = r(t) where r(z) is a rational approximation to the exponential 
e~ z given by: 


(1.27) r{z) = 1 - zb T {I + zA)~ 1 e. 

Expanding this expression shows that r(z) is a vtb. order approximation to the exponential if and 
only if (1.26) holds. Next, with regard to stability, an IRKM is said to be A q- stable if: 

(1.28) I r (*)| <1 Vz > 0, 
and strongly Ao-stable if: 


(1.29) sup |r(z)| < 1 Vzo > 0. 

z>z<3 

The former is required of all IRKM’s considered here, but for cases (1.1. ii) and (l.l.iii), the latter 
is assumed. Note that the spectrum of A, <r(A) is related to the poles of r(z) and in addition to 
the above, it is assumed throughout this work that: 
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(1.30) <r(A) C {z € C : tStz > 0, z ^ 0}. 

Finally, for cases (1.1. ii) and (1.1. iii), optimal results seem to require the mild condition that: 


(1.31) h 2 < ck 

unless Sh C Hq. The approach used is preferred since other attempts have led to (1.31) regardless 
of the boundary behavior of functions in Sh- 

Returning to the temporal discretization of (1.23), let a q-stage IRKM of order v > 1 be given. 
Then, the constants are well-defined by: 

v—1 

(1.32) £a jm m' = (- l) l l\ejA l e 0<l<u-l , 1 <j<q (m'l^^o = 1) 

m=0 

since their computation involves the inversion of the v x v Vandermonde matrix {m , }o<t»,i<i'-i- 
Next, with n*k = t* , and t n = nk, for v — 1 < n < n* — 1, suppose the approximations 
{U^}^ 0 C Sh are given, where U™ ~ u m and u m = u(x, t m ). Now, with L 2 = [L 2 Y, define 
the extrapolation operators £ n : [L^Y —* L 2 to have components: 

I /— 1 

(1.33) £?f h = £ VUj;- m ) 

m= 0 

with appropriate modifications for cases (1.1. i) and (l.l.ii). Also, with S* = define 

Ch : S/» -*■ S* and P 0 : L 2 -* S h by: 


£h = diag {L h } and 

qxq 


P 0 = diag{F 0 }- 

qxq 


Finally, let U£ +1 ~ u n+1 be given by what is henceforth called the principal scheme: 


(1.34) 


U£ = eU£-kAChUZ + kAP 0 £ n fh 
, UJ+' = U£ - kb T C h UZ + kb T P 0 £ n f h 


where C/£ € S/, and £ n fh respectively, are well-defined provided [I + kACh] is invertible and the 
approximations are sufficiently accurate. Here, ACh for example, is understood in the sense of 
composition of operators defined on Sf,. Note that (1.23) and (1.34) are only partially patterned 
after (1.24) and (1.25), i. e., the stages U£ are not fully implicit, and extrapolation circumvents 
the solution of a nonlinear system of algebraic equations for each n. 

Since the extrapolation for the principal scheme uses previously computed approximations, a 
starting procedure is required to generate Hence, for i > 1, define */,• = min(i/, 1 ), and 

as with (1.32), let the constants be determined by: 


i/,~i 


(1.35) Y, a ?m( m ~ n = HefA l e 


m= 0 


0 < l < Vi - 1, 1 < j <q ((m - = 1). 
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Then, with 0 < n < min(» - 1, v - 2), and my = min(i, v — 1), 1 < j < t, rrij = n,j = i, suppose the 
approximations {U™’ 3 }q C S'/, are given, where U™' 3 ~ u m . Next, define the extrapolation 
operators £ n >* : {LiY' — ► 1*2 to have components: 


(i.36) efh = 


m=0 


Vi-1 

E 

T7i=n+1 


E “gl/(*.‘“.pr'.vor')+ e or' -1 ), 

0 < n < i/,- — 2 


v,-i 


E » = .• - 1 < » - j 


\ m=0 

with appropriate modifications for cases (1.1. i) and (l.l.ii). Finally, let U£ +1,t — u n+1 be given by 
what is henceforth called the starting scheme: 


(1.37) 


f C/r = eUjf’* - kAZ h U^ + kAP^fn 


{ UJ? 1 ’* = U? - kb T C h U^ + kb T PaE"'' f h . 

With regard to the initial data U^' 3 = £7®, j > 0, for case (l.l.i), the following is sufficient: 


(1.38) Ul = P 0 «°. 

However, for cases (l.l.ii) and (l.l.iii), it is required that: 


(1.39) 


U% = [I + kL h }- l P 0 [I + fcL] u°. 


The cases (l.l.i) - (l.l.iii) are now analyzed separately in the next three sections. 


2 Semilinearities Independent of Spatial Derivatives 

In this section, the principal scheme (1.34) is analyzed for the approximation of the solution 
to (1.1) in case (l.l.i), and (1.9) is established. That the stages are well-defined depends on the 
Lemma below. Its proof involves a spectral argument after A is transformed to Jordan form, and 
the details are provided by Karakashian [13]. First, define the product space norm: 

ii*ii = {Eim 2 >* 

t=l 

Lemma 2.1 Provided (1.80) holds, [I + kACh] is invertible and: 

(2.1) WikZhYlI + kAUVxW <e||X|| VXeSfc, 0 <9<1. 


Now, sufficiently accurate starting approximations are assumed given and for v~l < n < n* — 1, 
an error equation relating Ufi — u> n to U £ +1 — o; n+1 appears below. For the sequel, make the 
definitions: 
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rj n = u n — w n , 


r = /(x,t",u"), 


«” = j2A l ed l t urk\ £?f = f£ = /(x, t", Ctf ), 

J =0 {=0 

£ = diag{L}, P B = diag{P £ }, r h = I - kb T £ h [I + k AC h ]~ l e. 

qxq qxq 

After some straightforward calculations, the following error equation is obtained: 


£ n+1 = r h £ n + kb T C h [I + kAC h ]- 1 P 0 {u n - eu" + JfcA£u" - kA£ n f } 


- kb T Z h \I + kAZ h \- x (P 0 ~ P E ){u n ~ eu") 


( 2 . 2 ) 


+ (Po - P E )(u n+1 - u") 

- P 0 {u n+1 - u n + kb T £u n - kb T £ "/} 

+ kb T [I + kA£ h ]- 1 P 0 [£ n f h -£ n f} = rh £ n + ^? + <f> n , 

l=i 


v- 1 < n < n* - 1. 

Stability follows with a spectral argument, and additional details connected with the following 
proposition are provided by Crouzeix [8]. 

Proposition 2.1 If the rational function ( 1.27 ) satisfies (1.28), then there exists a constant c < 0, 
such that: 


( 2 - 3 ) IkfcXll < (1 + cfc)||x|| Vx e S h . 

Furthermore, c < 0 if (1.29) is satisfied, and k is small enough. m 

Next, the order of consistency is established in the following. Also, <f> n is majorized by terms 
which are summed out in the convergence proof. 

Proposition 2.2 The terms {V'"}?=i °f (&•&) satisfy: 

(2.4) < c*(A ' + + IW«lk~ + lia,«l|r.„>. 

j=i 

Also, if o are given satisfying: 

( 2 . 5 ) 

then <j> n of (2.2) satisfies: 

(2.6) ||**|| < ctA'IWU + <*£ ll«~“l|. 

m= 0 


9 



Proof. The terms of (2.2) are considered in the order in which they appear. Since for 0 < / < i/, 
d\u eH 2 n Hi by (1.1): 

V V 

u n -eu n + kACu n = '£,A l ed\u n k l + ^A l+1 e[ L d\u n )k l+l 


l~l 


1=0 


1=0 


= — A v+1 ed% +1 u n k v+1 + k/i£,A l ed\j n k l . 

By (1.32): 

(2.7) 

where E has components: 

1 v - 1 

E I = E “*» / - *r '</(». “(»))<*»• 

Now with (2.1) and (1.1), it follows that: 

( 2 . 8 ) II^II < + ||3!7llo,„> < «v +1 {||ar‘«llo,o. + K»lk~>- 

Next, as with (1.32), the constants are well-defined by: 


1 / — 1 

£ n f = Y,A l t d \f n k l + E 
1=0 


^2 Pj™m l = l\v l ejA l e( 1 - 5 J0 ) 
m=0 


0 <l<v, 1 <j<q (m'l^^o = 1) 


where % is the Kronecker delta. Hence, define the extrapolation operator X n to have components: 

r;, 5 £Mi'+7) i<j<? 

m= 0 

so that for 0 < p < v, with ill-defined sums understood to be zero: 

X n u = Y,A l ed l t u n k l + F p 

1=1 


i ^ rt n l~ l 

^ £ Pim J tn V (t n + ^ ~ syd p+1 u(*)ds 1 < 3 < 7 - 


where F p has components: 

1 v 

E 

m=0 

Now note that is given by: 

t/>? = kb T [I+kACh]~ 1 PoC[u n -eu n - X n u] 

~kb T C h [I + kAC h }~ 1 {Po[u n ~ eu n - X n u] + [Po ~ Pe\X n u}. 
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(2.9) 


With (2.1), (1.19), and (1.20), it follows that: 

Mil < ck\\L[A' , ed v t u n k v - F''~ 1 ]\\ + c||P*'|| + c||[/> 0 - Pe\F°\\ 

< ck{k v + h r ){||d{'u|| 2 ,oo + ||ar 1 «|| 0 ,oo + || 9 ««||r,oo}. 

Now since: 

V »3 = / [Po ~ Pf;)d,u(s)ds, 

Jt n 

with (1.19) and (1.20), it follows that: 

(2.10) Mil < cA:/i r ||5tu|| ri00 . 

Next, using (1.1), (2.7), and (1.26), the following is obtained: 


v-l 


k l+1 


i/-i 


M = -Po{<? + I>i +lun 7nrTTT + 5> T ^[^l« n ]*' +1 - J2 bTAled tf nkl+1 - kbTE > 

{=0 \ l+1 )‘ 1=0 1=0 

= -Po{G + b T A v eLd(u n k v+1 - kb T E} 


where: 


Then using (1.1): 


1 / , < n+1 

G = — I (t n+1 - s)‘ / a^ +1 u(s)ds. 
Vi Jt n 


mi < ck^ipr- i «ii 0 ,oc + Hank*, + pr/iio,oo> 


(2.11) 

< cjfc*' +i {ii5r +1 «iio,oo+iiar«ii2 1 oo}. 

Now, (2.4) follows after combining (2.8) - (2.11). Finally, because of (2.5), (1.2), and (2.1): 


i/-i 

E 

trv^O 


v-l 


ini < ok 52 li/r m - r m n < <*,* y, nr— - »”- m n 


m= 0 


and (2.6) follows with (1.19). 

Now set: 

*(«) = HU + IMU + E liaHu + ll«r +1 «llo.» 

1= 1 

and (1.9) is finally established in the following. 

Theorem 2.1 Assume that (1.28) and (1.80) hold. Suppose are given satisfying: 


( 2 . 12 ) 


max 1 ||Dr-«"ll^e(* r + *■')>'(“)• 


Then , provided h N / 2 (h r + k l ') + 7o(h) is sufficiently small, {£/£}£_„ are well-defined by (1-84) and 
the following holds: 


(2.13) 


max \\U£ - u"|| < c{h r + k v )M(u). 

0<n<n* 
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Proof. Set 9{h,k) = h~ N / 2 (h r + k u ) + 70(h). It is first established that for 9(h,k) small enough: 

(2.14) .\\U? ~ W*\\ Loo < p. 

n 

Note that by (1.14), (2.12), and (1.21), for 9(h,k) small enough: 

HU? - « m ||ico < ch-^/ 2 ||ril + ||i? m || ico < C9(h,k) <p 0<m<u-l. 

Now suppose that for each h and k, there exists ann = n(h,k) such that v — 1 < n < n* — 1 and: 


(2.15) 

max ||^ - u'llr^ < p 

while: 


(2.16) 

\K +l -u n+1 \\ Loo >p. 


Given (2.15), inequalities (2.3), (2.4), and (2.6) can be combined for the error equation (2.2) to 
obtain: 


||£ ,+1 || < (1 + Ik) II^H + Cl k[h r + k v )M{u) + c 2 k J2 ||e'- m || u-l<l<n. 

m=0 

After summing this over v — 1 < l < n and applying the discrete Gronwall Lemma to the result, it 
follows that: 

(2.17) ||r +1 || < (1 + c„»)“~" +1 {ll£''~ l ll + + c„k £ ||ril> 

m= 0 

for cz = c -j- c 2 is > 0. Also, the exponential dependence on t* can be eliminated if cz < 0 [14]. Now 
by (1.14), (2.17), (2.12), and (1.21), for 6{h,k) small enough: 

II Etf +1 - «" +1 ||r 00 < ch- JV / 2 ||r +1 !l + ||»? n+1 ||Lco < c9{h,k ) < p. 

This contradicts (2.16), and hence, (2.14) is established. In fact, (2.13) follows from (2.12) and 
(2.17) after using (1.19). ■ 

See [14] for a proof of the following. 

Theorem 2.2 Assume that (1.28) and (1.80) hold. Then provided h~ N / 2 (h r + k 2 ) + 70(h) is 
sufficiently small, ore well-defined by (1.87) and (1.88), and the following holds: 

max. || Uf* - u> n || < c(h r + k i+1 )ff(u) 1 < * < v - 1. 

Therefore, (2.12) follows with U£ = U£’ v ~ 1 , 0 < n < v - 1. ■ 
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3 Semilinearities Depending Linearly on Spatial Derivatives 

In this section, the principal scheme (1.34) is analyzed for the approximation of the solution to 
(1.1) in case (1.1. ii), and (1.9) is established. Again, sufficiently accurate starting approximations 
are assumed given. Also, an error equation appears below, which differs in certain ways from (2.2). 
For the sequel, define the new terms: 

u" = Y,A l ed\v. n k l , g" = g(x, t n , u"), /" = g" • V«", g? = g(x, t", U£). 

1=0 

Note that here, u n contains one less term than its counterpart in (2.2). After some straightforward 
calculations, the following error equation is obtained: 


i/-i 


( n+1 = r h £ n + kb T [I + kAC h ]~ 1 PoC{u n - eu n - kA^2A l ed l t +1 u n k 1 } 


(3.1) 


1=0 


I /-1 


+ kb T [I + kAC h \-'Po{E "/ - '£,A l ed l t f n k '} 


1=0 


I/-1 


+ kb T [I + kAC^iPo - P E )^,A l ed l t +1 u n k l 


1=0 


1/-1 


{u n+1 -u n - kb T ^2A l ed l t +1 oj n k 1 } 


1=0 


+ kb T [i+kAc h ]- i p 0 {£ n f h -£ n f} = r*r+I>r+<r, 

i=i 

v — l<n<n* — 1. 

Now for the present section, stability is established in the norm: 

lllxlllis{(x,x)+*(£»x.x)}“ xes„. 

Proposition 3.1 Suppose that (1.29) is satisfied. Then for k small enough, there exists a constant 
c < 0 such that the following holds for ( S.l ): 


(3.2) \\\C + % < (1 + ck) mrili; + c*- 1 i:illV’rill? + c{\\(kL h rb n \\ 2 + \\(kL h )b n \\ 2 }- 

1=1 


Proof. Since it is assumed that v > 1, by (1.26) and (1.27), r(0) = 1 — — r'(0). So, let z\, 0 > 0 be 
chosen so that: 

\r(z)\ < 1 — Oz 0 < z < z\. 
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Next, define: 


_ 1 + (1 + e)z 
1 + z 

and fix e > 0 small enough so that with (1.29): 

\r € (z)r(z)\ < (1 + ez)(l — Oz) <1 0 < z < z\, 

and: 

|r £ (z)r(z)| < (1 + e) sup |r(*)| <1 zi < z. 

2>2l 

Then just as with (2.3), there is a constant c < 0 such that: 

(3.3) ||r £ (*L fc )r*xll 2 < (1 + cfc)|| X || 2 V X eS„. 

Now set c £ = 1 + e, operate on both sides of (3.1) with [I + c e kLh\ and integrate against f n+1 . 
Estimation of the resulting terms follows. 

With X n = [f + kLh]5$ n , it follows from (3.3) that: 

([I+c t kL h ]r k C, C +1 ) = (u(kL h )r hX n ,X n+1 ) < J(1 + S^IIIHg + 

By (1.18): 

([I + c e kL h ]W,t n+1 ) = ({k-h + c ( {kL h )^WM{I+LlH n+1 ) 

< cA:- 1 |||V’rlll?+^ll4e n+1 || 2 l</<4. 

Finally: 

((/ + <:««»]*", f" +1 ) = ({(t£»)-i+c.(*I»)*W",(*£»)ir fl ) 

< c{ii(H, k )-4"ii“+ii(tL ft )i]^r}+Hi4{" +i ii 2 . 

Combining the above inequalities: 

(r +1 , e n+1 ) + (i + e)k(L h c + \ c +1 ) < |(i + skiing + liiir +1 iu 2 i + ^-‘Dwiig 

i=i 

+c{\\(kL h )-h n \\ 2 + ||(^) 4 "|| 2 } + ek\\LU n+1 \\ 2 

and (3.2) follows readily. ■ 

Now an analogue to Proposition 2.2 must be established for the stronger norm. 

Proposition 3.2 Assume that (1.81) holds. Then the terms {ip”}f =1 of (8.1) satisfy: 

(3.4) e \\m\\i < ck(h r + fc ,/ ){ii 5 r l «iio > oo + iiar«n 2j oo + phu}. 

j=i 
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Also, if {U£ m }m=o are given satisfying: 


max \\U?~ m - u n - m \\ Loo < p, 


; max_i||?7” m \\ w i^ < c. 


for some c > c u = sup ||u(t)||jyi,oo, then <f> n of (S.l) satisfies: 

o<t<r 


(3.7) ||(*U)-4“I| ! + < c*(h’-|l“llr.«o) 1 + ek~£ llle n - m |lll- 

m= 0 

Proof. By (2.1): 

(3.8) IH^IIl! < ck\\tA v ed v t u n k v \\ < cfc I ' +1 ||ar«|| 2l oo 
Next, by (1.32), an analogue to (2.7) follows. Hence, by (2.1) and (1.1): 

( 3 . 9 ) in*? it, < ^ ,,+1 iiar/iio, t » < »*- +i {iiar +1 «iio.» + U8?«ik~}- 

Now, recalling the development prior to (2.9), with (2.1) and (1.30), the following is obtained: 
Kill, < c||(^ - Po)[A-ed : 'jW - /- 1 ]|| + c|| (P E - fl,)F°||. 

Using s — 2 and s = r in (1.19) and (1.20), with (1.31) it follows that: 

(3.10) Hl^llli < + ckh'\\d t u \ | r)0O < ck{k v + h r ){Ku|koo + ||d««||r,oo} 

Next, set: 

H = u n+1 - u n - kb T A l ed\ +1 u n k l 

1=0 

so that t />4 = PeH and: 

Kill? < ll«ll ! + *IIW?IKII £ + |t ! ||Po£fl'll s 

< c\\Hf+4(P c -I)H\\*+ck*l\H\\l 

This estimation is then continued using s = 2 in (1.19) and applying (1.31): 

lltf 111 < c||^ II + c(h 2 + k)\\H\\ 2 < c\\H || + <*||ff||,. 

By (1.26): 

H = G p + Y2 d\u n j - J2 b T A l ed l t +1 u n k l+1 = G p - {u - p)d^u n ^ p = v - 1, v 


where: 


1 r* n+1 

G P = 1 / (t n+1 - s) p d p+1 u{s)ds. 

P ' Jt n 
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Hence: 


(an) imiii. < ^lie'll + <*iio*“ I - ar«"£ii> < -* ,,+1 {iiar +1 -‘iic..~ + ii«r»ii2,»>- 

Now, (3.4) follows after combining (3.8) - (3.11). Turning then to the Gronwall terms, by (1.17) 
and (2.1): 


(3.12) 


m=0 
^L 1 1 

+ C ^E ||7^g n-m • V(i7^* _m — u n_m )|| 2 


m= 0 

From (3.5), (1.3), (3.6), and (1.19), it follows that: 


0 = ±i 


ii(gr w -g n_m )-vc7r m ii < cc,iie n - m -i? n -"‘iiiicfr m ik».- 

(3.13) 

< cc||f n m || + ch r ||u|| r>O0 0 < m < v - 1. 

Next, set: 

D n-m = g n-m . V(f/” _m - u n ~ m ) 0 < m < U - 1. 

By (1.8) and (1.12) with s = 2: 

\\T^ D n ~ m \\ 2 = \\T*D n - m \\ 2 + {[T h -T}D n - m ,D n - m ) < c||£ n - m ||ii + ch 2 \\D n ~ m \\ 2 . 

Using (1.6) and (1.19): 


|Z? n-m ||_ 1 = sup 


|([ur m - u"-"*],V . [<pg n ~ m })\ 

Mi 


< cc g \\r~ m - n n ~ 


< c|K"-’*|| + cfc r ||tt n ^ m || r . 

By (1.18), (1.31), and (1.19): 

< c Ctf h|ir- m - r? n - m iii < c|ir- m iiii + ch r i|u”— i| r . 

Combining the above inequalities: 

(3.14) HT/g"-” 1 • V(ur m - U n_m )|| < c ||| r - m ||| 1 + ch'\\u\\ r>00 0<m<u-l 

and (3.7) follows from (3.12) - (3.14). ■ 

Following (1.31), it is claimed that the relation between h and k can be avoided if Sh C Hq. 

This can now be seen in the preceeding proof. If the range £(T fc 2 ) C Hq, then using (1.17), an 

1 

inequality such as (1.8) can be established for T £ . Given the latter, it would not be necessary 
to triangulate with T in order to obtain (3.14). Now if this approach is taken and (1.31) is not 
assumed, then the consistency result is: 

4 

EllWllll < ck(h r + k v + h 2 k v ~ X ){ K +1 tx|| 0 .oo + liaHkoo + IMIr.oo}- 

1=1 
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However, since (1.31) is such a mild condition compared to the requirement that C Hq, the 
details of the result suggested here are not provided. Instead, (1.9) is established as follows. 

Theorem 3.1 .Assume that (1.29), (1.S0), and (1. SI) hold. Suppose {U™}'£l) 0 are given satisfying: 


(3.15) 


maxjDT - " ro Hli < W + k v )fl { «). 


Then provided h N ! 2 k 1 / 2 (h r +k v ) + 71 (h) is sufficiently small, {{/£}"=„ are well-defined by (I.S 4 ) 
and the following holds: 


(3.16) 


o maxJ|^-«l<c(h- + fc‘') J V(«). 


Proof. Set 0o(h,k) = h N / 2 {h r +k' , )+io(h), Oi(h, k) = h N ! 2 k 1 / 2 (h r + fc*')+ 7 i(h), and let c > c u . 
It is first established that for 9o(h,k) and 9i(h,k) small enough: 


(3.17) 


max 110 ? - « 1 ioo < p 


0<m<r»* 


max WUJTWwi,- < c. 


0<m<n l * 


Note that by (1.14), (3.15), and (1.21), for 9o(h,k) small enough: 

ns? - « m llL_ < or^nril + IIVIIi. < «*o (*,*) < P 

Also, by (1.15), (1.18), (3.15), and (1.22), for 9i(h,k) small enough: 
WUTWwu- — + ll , 7 m ||» rl - 00 +c u < c9\(h,k) + c u < c 


0 < m < v — 1. 


0 < m < u — 1 . 


Now suppose that for each h and k, there exists an n = n{h,k) such that v — 1 < n < n* — 1 and: 


(3.18) 

while: 


max \\ U h “ < P> 


0 <J<n" 


max 1117^11^1,00 < c 


(3.19) either \\U^ +1 - u n+1 \\ Lca > p, or > c. 

Given (3.18), inequalities (3.2), (3.4), and (3.7) can be combined for the error equation (3.1) to 
obtain: 


ie ,+1 m? < ( 1 + 5*)k'iiiS + «!*[(*' + nswi* + «»t E ik'-h; 


v — 1 < / < n. 


After summing this over v — 1 < / < n and applying the discrete Gronwall Lemma to the result, it 
follows that: 

( 3 . 20 ) iur K nii < (i + c 3 fc) n - ,/+1 {|iir _1 iii? + cit*[(h r + k^Miu)] 2 + c s k £ inmg} 

m= 0 

for C 3 = c + c%v > 0. Also, the exponential dependence on t* can be eliminated if C 3 < 0 [14]. Now 
by (1.14), ( 3 . 20 ), (3.15), and (1.21), for 9 0 (h,k) small enough: 

||t/; +1 - < c(.- K / ! ||r +1 ll + lli“ +, IU. < cHo (h,k) < p. 
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Also, by (1.15), (1.18), (3.20), (3.15), and (1.22), for 0i(/i,Jfc) small enough: 

ll^h +1 ||w l -°° < c h A ^ 2 ||f n+1 ||i + ||'7 n+1 ||w ,1 .°° + c u < c0i(h,fc) + c u <c. 

This contradicts (3.19), and hence, (3.17) is established. In fact, (3.16) follows from (3.15) and 
(3.20) after using (1.19). ■ 

See [14] for a proof of the following. 

Theorem 3.2 Assume that (1.29), (1.S0), and (1.81) hold. Then ifh~ N l 2 k~ l l*(h r + A: 3 / 2 ) + 71 (h) 
is sufficiently small, are well-defined by (1.87) and (1.89), and the following 

holds: ~~ 

n< <™ a r ,1 “ ""Ml - C ^ r + *‘ /2+1 )^( u ) 1 < * < 21/ - 1. 

Therefore, (8.15) follows with Uj) = 0 < n < v — 1. ■ 

4 General Semilinearities 

In this section, the principal scheme (1.34) is analyzed for the approximation of the solution to 
(1.1) in case (1.1. iii), and (1-10) and (1.9) are established. For this, the basic structure of previous 
sections is followed and in fact, several of the estimates of section 3 are readily adapted. First, the 
error equation (3.1) serves here but with: 

r = /(x,t",u n ,Vu n ), fZ = f( X ,t n ,UZ,VUZ). 

Stability is now established in the norm: 

lllxlll 2 = {{ L hx,x) + H L lx,x)}* x^s h . 

Proposition 4,1 Suppose that (1.28) is satisfied. Then , for k small enough, there is a constant 
c < 0 such that: 

( 41 ) mr +1 iii2 < (i+2fc)iiiriii2+cfc- 1 ii[/+fcL fc ](r +1 -r fc e n )ii 2 . 

Proof. With c € = 1 + e and x n = [Lh + f rom (3-3) it follows that: 

([£fc + c,kLl\r h (", £" +1 ) = (r,(H») r/ l x n ,x” +1 ) 

< |(i+«)inriii? + |iiir +, iiil 

Also: 

([Z, + c ( kLl\(C +1 - r h C), r +1 ) < cA:- 1 ||[J + kL h ](C +1 ~ rfc £ n )|| 2 + ^Ill £ n+1 ||| 2 

and (4.1) follows after summing these inequalities. ■ 

Now, the new terms of (4.1) must be treated differently. 
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Proposition 4.2 Assume that (1-Sl) holds. Then the terms {V*p}?=i satisfy: 


(4.2) 0(/+*L»l*r|| < c k(h r + fc"){||3)' +1 u||o J oo + ||ar«lll.~ + IMU)- 


1=1 
Tn — m \ i /— 1 


Also , if {U£ m }JUo ar€ 0 iven satistying: 
(4.3) 

then <f> n satisfies: 


max II U? m — u" m ||^i,oo < p , 

0<m<«/— l" h " W 


(4.4) 


III/ + *£»]« < cfcfe r - 1 ||«||r,oo + ||r-"ll,. 

m— 0 


Proof. For (4.2), inequalities (3.8), (3.9), and (3.10) are readily extended by using 9 = 1 in (2.1). 
Then, with = PeH as in Proposition 3.2: 

||[/+ kL k M\\ < ||V*4 II + k\\P 0 LH\\ < || (P E ~ I)H\\ + \\H\\ + k\\LH\\ 
and the remaining component of (4.2) follows as with (3.11). Now, by (2.1): 

u(«„)V"ii < ck'£]\nuz-’',vurn - /(«-“. vqr"ii 

m= 0 


+ - /(«”-”. vp- 

m=0j=l 


e = o,i 


where V 0 n = VU£, Vj} = Vu n and: 

Vp={d Xl d x .u n , d Xj+1 U^...,d XN UP) T l<j<N-l. 

Using (4.3), (1.4), (1.19), and (1.18): 

u/(tr m ,vur m ) - /(« n - m , vur m )ii < cc,iie n ~ m - f? n - m n 

< c|||e n - m ||l 2 + c/i r ||«|| ri oo 0 < m < v — 1. 

Similarly: 

|| / ( u n-m,^n_-m)_ / ( u n-m,^»-m)|| < CC p 1 1 ^ - fj n ~ m )ll 

< c|||e n - TO ||l 2 +^ r “ 1 |l«llr.oo 

1 < j < N, 0 < m < v — 1. 


and (3.7) follows after combining the last three inequalities. ■ 

The groundwork for an H 1 estimate is now complete. For an L 2 estimate, the natural impulse 
is to press the details surrounding (3.14) for a generalization to the case (1.1. iii). In search of an 
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analogue to the H 1 estimate prior to (3.14), it is tempting to suppose the existence of constants 
p> c p > 0 such that for 0 < t < t*, the following are satisfied: 


and: 


r W U U 2 e W 1 ’ 00 satisfying max \\d x .(U m - u)|| ioo < p 

| tn — 

||/(t,u(t), VU 2 ) - VI7 1 )||_ 1 < c p \\U 2 - lh || 


{ Vv G W 2 '°° satisfying ~ «)l| < P 


f(x,t,u(x,t),Vv(x)) G W 1 ’ 00 . 

However, it is shown in [15] that these conditions are actually equivalent to the following: 

3 {/ m (x,t,u(x,t))}£ =0 c W 1 *, f=(fuf 2 ,...,f N ) T 

such that Vtx? G W lt °° satisfying max^d^. (w - u)!!^ < p 

f(x,t,u{x,t),'Vw(x)) = fo(x,t,u(x,t)) + f(x,t,u(x,t)) • Vu/(x). 


In fact, this is at the heart of the division (1.1. i) - (1.1. iii). Nevertheless, it is clear that a new 
approach is required for the term <f> n before an optimal L 2 estimate can be proved. 

Proposition 4.3 Assume that (1.81) holds and suppose 0 are given satisfying: 


(4.5) 


max ||{7? m - u n m \\ W i, ao < p, 


Then with e G (0, 1), ppf >e = maa:{0, y — 1 + e}, <j> n satisfies: 


v-l 


||(ti k )-4”l| J + ll(Mj4 n || 2 < c*('>'|MU) ! +<*£|K n -’”lll? 


(4.6) 


* 71=0 


V— 1 

+ c(e)kh~^ N < € ^2 \\U£~ m — u n_m ||i. 

m=0 


Proof. By (1.17) and (2.1): 


i/-i 


\\(kL h yr \\ 2 < ck^\\f{U£- m ,VUZ- m )- f{u n - m ,VUZ- m )\\ 2 


(4.7) 


m= 0 

V-l J 

+ ck J2\\Th[f(v n ~ m ,VU£- m )- f{u n - m ,Vu n ~ m )}\\ 2 e = ±i. 

m= 0 


From (4.5), (1.4), and (1.19), it follows that: 


20 


(4.8) 


Muz- m ,vuz~ m )-/(«"-”*, v^- m )ll < cc /) ||e n -” l -»/ n - m || 

< c||e n - m || + c^||u|| r>c 


0 < m < v — 1. 


Next, combining (4.5) and (1.5), ggjT 1 and are well-defined by the 

following: 




/(u n " m , VU£- m ) ~ f{u n ~ m , Vu n-m ) = E^W””* - u n ~ tn )dd j f(u n ~ m , Vu n-m )+ 

y=i 

E 5 It .(^r m - « n ~ m )d Xj (ur m - u n - m ) f\i-s)di dj n« n - m , vw n - m + °w~ m - « n_m )]^ 


*.i=i 


(4.9) 


W AT 

= E-°r m + E ^,” rm 0 < m < */ - 1. 

i=i l 


Now, by (1.8) and (1.12) with s = 2: 

l | T ; fi >?- m || 2 < || T 5£>;- m || 2 + ([ T fc - £>?- m ) 

< c||Z?;- m ||?. 1 + ch 2 ||£>?- m || 2 l<j<N, 0<m<u-l. 
The first of these two terms is estimated as follows using (1.7), and (1.19): 

l|D > "- 1 - S iMil 

< CCf\\r- m -V n - m \\ < C|||r-” l ||| 1 + ch’-||ti ,l - m ||r 

1 < j < N, 0 < m < v - 1. 

For the second term above, with (1.7), (1.18), (1.31), and (1.19) it follows that: 

h\\D^- m \\ < cc / h||e"- m -r,"-” l || 1 < c|||e- m ||| 1 + ^||u'*-~|| r 

1 < j < N, 0 < m < v — 1. 

Combining the last three inequalities: 

(4.10) ll^i>r m |l ^ c lll^" m |lli + ^IMkoo 0 < m < 1 / — 1. 

i = i 
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Next, using (4.5), and (1.5): 


(4.11) 


i \(E?r m T?<p)\ 

iit^-hi = sup « ; 


(pGL 2 


< CC„\\UZ m - U n m ||i sup 


v>eL 2 Ill’ll 


l<i,j<N, 0 < m < i/ — 1. 

By the Sobolev Imbedding Theorem (See Adams [1]), since Sj, C W 1,0 °, for e e (0, 1): 


HTMioo = sup |T^Xx)| < c(6)||T fc V|| V yl,W(l- < )-‘ € L 2 . 

xsn 


In case iV = 1, with e = | this is sufficient for what follows. For N > 2, note that by (1.16): 

< ch- N /> +1 -‘\\T*<p\\i V<p E L 2 . 

In any case, from (1.17) and the last two inequalities, it follows that: 

(4.12) II^IUoo < c(e)h-^|M| V<p € Li. 

Combining (4.11) and (4.12): 


(4.13) \\Tji J2 E ?f m \\ < c{e)h-P^\\U^~ m - u— 1|? 0 < m < v - 1. 

U=i 

Then, combining (4.10) and (4.13) for (4.9): 

rJ[/(«~“,VV*-“) - V«“)]|| < eVWU + clle—11, 

( 414 ) + c(()h~n<ui7"- m - 


0 < m < v — 1. 

Now, (4.6) follows after combining (4.14) and (4.8) for (4.7). ■ 

Note that the remarks made prior to Theorem 3.1 apply here as well. Nevertheless, the conver- 
gence results (1.10) and (1.9) are now established as follows. 

Theorem 4.1 Assume that (1.29), (1.S0), and (1.81) hold. Suppose are given satisfying: 

( 4 .1 5 ) - <,”- 111 , < + rwu). 

Then provided h~ N ! 2 {h r ~ 1 + k v ) + 7i(h) is sufficiently small, {{/£}£!,, are well-defined by (1-84) 
and the following holds: 

(4.16) o max J|C£ - u*\\x < c(h^ + k')M(u). 

If also satisfy: 
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(4.17) 
then: 

(4.18) 




n>M.|TC-«l<4(A' + t-)41(,.). 


0<n<n’ 


0<T< X n* l|C/r ~ ^ 


Proof. Set 0(h, k) = h N ! i (h r 1 + k u ) + 7 i(A). It is first established that for 0(h, k) small enough: 
(4.19) 

Note that by (1.15), (1.18), (4.15), and (122), for 0(h,k) small enough: 

\\UJT - « m |Ui,oo < ch-^' 2 ||e m ||i + !|f? m !! W ri.oo < C 0 (h.k) <p 0 <m<u-l. 

Now suppose that for each h and k, there exists an n = n(h,k) such that v — 1 < n < n* — 1 and: 

l ..l i 


(4.20) 

while: 


max IIP* - « \\w'*° < P » 

USfSn 


(4.21) \\Up 1 -u n+1 \\ wl , 00 >p. 

Given (4.20), inequalities (4.1), (4.2), and (4.4) can be combined to obtain: 

i/-i 


ine‘ +, ig < (i + 3*)iii€'iiil + «iM (*'-■ + D mf 


'-“US 


- 1 < / < n. 


m=0 


After summing this over u — 1 < / < n and applying the discrete Gronwall Lemma to the result, it 
follow that: 


i/-i 


(4.22) Hir'-'III! <(! + «*) 


n-»-+l,|||ei/-l||Lt 


„(•(((.'-■ + + csk £ iiiriiii) 


m=0 


for C3 = c+ C 2 V > 0. Also, the exponential dependence on t* can be eliminated if C3 < 0 [14]. Now, 
by (1-15), (1.18), (4.22), (4.15), and (1.22), for 0(h,k) small enough: 

U^n+l _ u n+l| (Hrl oo < c h- JV / 2 ||r +1 !!l + ||»7 n+1 !kx.« < C$(h,k ) < p. 

Hence, (4.21) is contradicted and (4.19) is established. In fact, (4.16) follows from (4.22) and (4.15) 
after using (1.19). 

Now, if (4.17) holds, (4.18) is obtained as follows. Given (4.19), inequalities (3.2), (4.2), and 
(4.6) can be combined to obtain: 


\H l+1 W\l < (1 + ck)\\\Z% + c 4 k[{h r + k^Miu)} 2 + c 5 k J2 llle'-HII 2 ! 

m= 0 


*/ - 1 < / < n* - 1 
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since for small enough e and $(h, k), r — 2 — p# |£ > 0 and: 


(c(h r + k v )M{u) + c{e)h-P^[{h r ~ l + k^Miu)} 2 } 2 < c[h r - 2 ~PN.‘ + 0{h, Jfc)] 2 [(h r + k v )M{u)} 2 

< CiKh' + k^Miu)} 2 . 

After summing over v — l</<n<n*-l and applying the discrete Gronwall Lemma to the 
result, it follows that: 

(4.23) nir +I iii! < (i + «6i)”-' +1 {ine‘" 1 in? + «<»•[('*' + nm ? + £ iirng> 

m= 0 

for c$ = c + CiV > 0. Again, the exponential dependence on t* can be eliminated if eg < 0 [14]. 
Finally, (4.18) follows from (4.23) and (4.17) after using (1.19). ■ 

See [14] for a proof of the next Theorem. It requires the inequality: 

(4.24) Il4x|| < efc^Hxll V X eS„ 

which depends on inverse properties such as (1.14) - (1.16). 

Theorem 4.2 Assume that (1.29), (1.S0), (1-81) and ( 4 - 24 ) are satisfied. Then provided 
h~ N ! 2 {h r ~ 1 + k z ! 2 ) + 71(h) is sufficiently small, {^h’*}o<lKmin(t v-i) are well-defined by (1-87) 
and (1.89) and the following holds: 

max IlltC 4 - « n ||| a < c^" 1 + Jfc </2+1 )JV(t») 1 < i < 2u - 1. 

Therefore, (4-18) follows with Ufi = 0 < n < v — 1. Furthermore: 

max 1 1| 17?’* - u n \L < c(h T + k^^Jdfu) 1 < i < 2v - 1 

and (4.17) follows with U£ = U^’ 2v ~ l , 0 < n < v - 1. 

5 Computational Aspects. 

Note from (1.34) and (1.37), the need of an efficient method for solving linear equations of the 
form: 


[I+kACh}* = $ V^eSh. 

In connection with parallel implementation, it can be seen below that the preferred methods are 
those for which the eigenvalues of A are distinct. So, let these be called multiply implicit (MIRK) 
methods. 

First, consider the use of IRKM’s which have high order with respect to q, such as the Gauss- 
Legendre (v = 2 q) or the Radau [u = 2 q- 1) methods (See e. g., Dekker and Verwer [9]). These are 
Ao-stable, and the latter are strongly Ao-stable while the former are not. Also, though each class 
of methods satisfies (1.30), for neither is o-(A) C R. Nevertheless, these are (complex) MIRK’s, for 
which A can always be transformed to quasidiagonal form: 
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A = S' 1 AS 


where for some m, 1 < m < q: 
and: 


A = diag {A<}, 

l<t<m 


either A,- = Af > 0 or A,- = 
Then, the linear system can be written as: 


a,- -fit 

Pi a i 


cti, Pi > 0 1 < * < m. 


[I+kA£ h }(SV) = (S«) 
which decouples to equations of the form: 

[I + k\Lh]xp = <j> rp,<f>€Sh 


and: 


I + kccLh 

~kpL h 

V> i 


<t> i 

kPL h 

I + kaLh 

. ^ 2 . 


. 


V’i) 4*i £ Shi *' — 


Note that the subordinate equations can in principle, be solved simultaneously. Further, the 2x2 
system is equivalent to the following pair of equations: 


[I + kaL h ]rpi = </>i + kpL h rp 2 


and: 


[I + 2 akL h + (a 2 + P 2 )k 2 L 2 h ]tp 2 = [I + kaL h \<j> 2 - kpL h <j> x = xo- 
Following Baker, Bramble and Thomee [2], since for x E R: 

Zi = 2a + P~ l {a 2 + P 2 )i, z 2 = a + pi, 

complex arithmetic can be used to obtain: 


1 + 2ax + (a 2 + P 2 )x‘ 


= 1 


Z IX 

U + z 2 x . 


V "2 = Xo - + z 2 kl] x xo}. 


Next, in spite of the order barrier v < q-\- 1, consider methods for which A is similar to a matrix 
of the form: 

Ai 0 , 

A,- > 0, !<»<? 


A = 


0 2 \ 2 
0 


0q A, 


$i = 0 or 1, 2 < i < q 
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so that the block system reduces to equations which are linear in L^: 


[/ + kX\Lh] V*i — <f> 1 

[I + kXiL h ] ipi = 4>i ~ k6iL h ipi- 1 , 2<i<q 


ipjyfyeSh, 1 <j<q. 


A class of (real) MIRK’s is obtained in case A,- 7 ^ Ay, t 7 ^ j and 0 t - = 0, 2 < » < q. The 
depth of decoupling these methods allow in the above equations makes them appear, in a parallel 
environment, more attractive than their complex counterparts above. However, the cost of this 
advantage is reduced order. 

On the other hand, the class for which A,- = A, 1 < t < q is said to contain singly implicit 
(SIRK) methods. Since these require the formation and factorization of only a single matrix with 
the dimension of S SIRK’s are preferred if computations are to be performed on a serial machine. 
A selection from the previous set was made for the example discussed below. 

The following problem is of the class considered in Section 2: 


d t u = d\u + f(x,t,u) in (- 1 , 1 ) x [ 0 , . 1 ] 

u = 0 on {- 1 , 1 } x [ 0 ,.l] 

u(x, 0 ) — 1 - x 2 in (— 1 , 1 ) 


where: 


f(x,t, u ) 


2 - (1 - x 2 ) log(l + X 2 ) 2t(l - 5x 2 ) 

(1 + x 2 y (1 + X 2 ) t+1 


4x 2 t(t + 1)(1 — x 2 ) 

(1 + x 2 y + 2 ’ 


X 2 = 1 - (1 + x 2 )‘u. 


The solution is given by: 


u(x,t) 


1 — X 2 

(1 Tx*Y' 


For the spatial discretization, the Ordinary Galerkin Method was used and S* was constructed of 
smooth cubic splines defined on a uniform mesh. For the temporal discretization, the following 
three-stage diagonally implicit (DIRK) method was used: 


7 

0 0 

7 


7 0 

1 

2 

2 7 

1-47 7 

1-7 

P 

1-2/3 fi 



'»= ?+ 75 c “(n) 


In the sequel, let (1.34), (1.37), and (1.38) be identified as the modified method. For comparison, 
a classical method is now introduced. Define T n : L 2 — * L 2 to have components: 


//*($) = /(x,f" + fc7-,-,<k) $eS fc . 


Next, let U%'° 
take: 


eZ/£ = ePou 0 and after the stages {U£ are computed as indicated below, 

min(n,i/) 


op = E (-i) m+1 " 0;- 

m= 1 V / 


1 < n < n* - 1. 
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k,h 

CPU Time (sec) 

Z/2 error (xlO 9 ) 

Order 

1/50 

12 

21.4 


1/60 

17 

9.94 

4.20 

1/70 

23 

5.30 

4.08 

1/80 

31 

3.09 

4.03 

1/90 

39 

1.93 

4.01 

1/100 

45 

1.26 

4.00 


Table 1: Modified method 


k,h 

CPU Time (sec) 

L 2 error (xlO 9 ) 

Order 

1/50 

51 

408. 


1/60 

67 

245. 

2.80 

1/70 

87 

158. 

2.84 

1/80 

112 

107. 

2.88 

1/90 

136 

76.4 

2.91 

1/100 

168 

56.0 

2.94 


Table 2: Classical method 


Then define: 


V? = [I + kAZ h \~ x {eUZ + kAPoT^Uf- 1 )} l > 1, 0 < n < n* - 1 

and: 

fjn = jj«,max(v + l-n,l) 0 < „ < „* - 1. 

Finally, take: 

U£ +l = U£- kb T Z h U£ + kb T P 0 T n (U£). 

Now, in addition to the modified and classical methods, let a hybrid method be given by (1.34), 
(1.37), and (1.38), but with l\A l e replaced by T*e in (1.32) and (1.35). 

These three methods were tested on the ICASE SUN 3/180. In Tables 1-3, the Lz errors: 

E(k) = E(k,k), E(h,k) = |W-«”'|| 


are reported together with estimates of the convergence order obtained according to the formula: 


Order = 


\og{E{k 2 )/E(k 1 )) 

log(fc 2 /fci) 


As expected, the classical method requires much more time because of the additional function 
evaluations. Also, for the present example, the classical and hybrid methods offer roughly the same 
accuracy. However, at least according to theoretical efforts, this is not to be generally expected. 
Finally, note that the literature apparently does not contain an explanation of the better than 
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k,h 

CPU Time (sec) 

Li error (xlO 9 ) 

Order 

1/50 

12 

402. 


1/60 

17 

244. 

2.75 

1/70 

22 

158. 

2.82 

1/80 

29 

107. 

2.87 

1/90 

37 

76.3 

2.91 

1/100 

46 

56.0 

2.94 


Table 3: Hybrid method 

second order convergence seen in Tables 2 and 3. In particular, with Lv = d 2 v, v € H 2 n Hi 
the above solution is not even in the domain of L 2 . Nevertheless, only second order convergence is 
demonstrated for example, in Experiment 7.5.2 of Dekker and Verwer [9], where a stiff nonlinear 
ordinary differential equation is considered. Further, the modified method has been applied to this 
problem to give not only fourth order convergence, but accuracy exceeding that reported for any 
method discussed in the Experiment. 
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